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Abstract 

We present a general procedure for calculating the partition function of an Ising Model 
on a one dimensional Fibonacci lattice in presence of magnetic field .This partition func- 
tion can be written as a sum of partition functions of usual open Ising chains in presence 
of magnetic field with coeffecients having Fibonacci symmetriesAn exact expression for 
the partition function of the usual open open Ising Model is found. We observe that "Mir- 
ror Symmetry" is a characteristic property of all Fibonacci generations . Further nth 
and (n + Q)th generations have the same topology. We have also established a recurrence 
relation among partition functions of different Fibonacci generations from nth to (n+6)th. 
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Introduction 

Studying Ising model on a Fibonacci chain in presence of magnetic field is interesting 
in its own right since no proper analytical procedure exists for evaluating the partition 
function. However scaling forms of thermodynamic functions for such system have been 
studied using renormalization group technique [1] through one step decimation. Here we 
reformulate the expression for the partition function by breaking the transfer matrix into 
two particular noncommuting matrices . This formulation enables us to calculate the 
partition function for the usual open Ising chain. The result is quite nontrivial in contrast 
to the expression for the partition function of the closed one [2]. The same formulation 
helps us to express the partition function of Ising model on Fibonacci chain in presence of 
magnetic field as a sum of partition functions of usual Ising open chains with coefficients 
containing Fibonacci symmetry. We have also studied some symmetry properties of the 
Fibonacci chain. Using a special symmetry property ("Mirror Symmetry") and the usual 
trace map relation we have established a recurrence relation among the partition functions 
of different Fibonacci generations. This includes all the partition functions starting from 
nth up to (n + Q)th generations. We observe that mirror symmetry is a characteristic 
property of each Fibonacci generation with nth and (n + 6)th generations having same 
topology . 



Exact partition function for open Ising chain with magnetic field 

The one dimensional Ising model consists of a chain of N spins Si = ±1 
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;i = 1,2, , N with nearest neighbour interactions e^+i. The 

Hamiltonian is given by: 

N-l N 

7~L — — £j,i+iSiS i+ i — Hy^Sj (1) 
1=1 1=1 

For a uniform lattice e iji+ i = e,the partition function is given by: 

Z° N (T, H)= £ f(Su S 2)f(S2, S s ) f(S N -i, S N )f (S N , S l ) (2) 

Si,S2,----,Sn=— 1 

with /($, S i+1 ) = exp\j3eS i S i+1 + \(3H{S, + S i+1 )\; f (S N , S } ) = [f(S N , 5i)] e=0 . Here 
the superscript o stands for the chain with open boundary condition. Therefore the parti- 
tion function (2) can be written in terms of transfer matrix as: 

Z° N (T,H) =TrP N - 1 P (3) 

where 

P = ^(l + V = ^ 1 (l + ^) (4) 

P = [P]e=0 = (1 + A) (7i = <7x(l + A T ) (5) 

with r = exp(-2(3e) , A = ( 6 Q ) , A T = ( J H 6 Q )andai=^ J ) . 

The general formula of the partition function for even and odd number 

of spins (i.e., odd and even number of bonds) can be derived by using equations 

(3) and (4) as : 



and 



where 



Z° 2N (T, H) = r N "Tr{\ + Xl )(l + x 2 )....(l + x 2N _ 1 ){\ + A T ) (6) 
Z° N+1 (T, H) = r N Tr(l + Xl )(l + x 2 ) (1 + x 2N )(l + X)a, (7) 



A A T 

%2i+i = -, %2i = — ;i = integer (8) 

The above equations show that A, A T are the signatures for the transfer matrices 
corresponding to bonds in odd and even positions. In the case of a chain with closed 
boundary condition the last factor in eqn.(2) is /(Sat, Si) and consequently the partition 
function takes the form: 

Z C 2N (T, H) = r N Tr{\ + zi)(l + x 2 ) (1 + x 2N ) (9) 

Z 2N+ i = r N+l >Tr(l + Xl )(l + x 2 ) (1 + x 2N+1 )a Y (10) 

Here the superscript c indicates closed chain. One can show by elementary calculation 
that eqns.(9) and (10) reduce to the well known form[3] 

Z C N {T,H) = \ N + +\ N _ (11) 

where 
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A± = r~ 3 [cosh (/3H) ± yf(swh 2 (PH) + r 2 )] (12) 

are the eigenvalues of the transfer matrix P. The expression for the partition function 
in the case of an open chain with even number of spins can be derived from eqn.(6) as 
follows: 



Z° N (T, H) = r N - l iTr(l + Xl )(l + x 2 ) (1 + x 2N _ 1 )(l + A T ) 

= V^Z C 2N (T, H) + r N ~Hl - r)Tr{\ + Xl )(l + x 2 ) (1 + x 2N _ 1 ) 

= V^Z C 2N (T, H) + ^(1 - r)Z c 2(N ^(T, H) + r N ^(l - r) 

xTr(l + Xl )(l + x 2 ) (1 + x 2N _ 2 )x 2N ^ (13) 

The last term in the above expression can be written in terms of the eigenvalues of 
the transfer matrix P viz.A±. By following the method of induction : 



r N -\(\ - r )Tr{\ + Xl )(l + x 2 ) (1 + x 2N -2)x 2N -i 

= (i _ r )r N -^cosh 2 ((3H) ]T (^-) 



i=0 



\ 2(/V-l) _ x 2(iV-l) 

4(1 - r)r^ cosh 2 ((3H) + 2 — ^ ( 14 ) 



So eqn.(13) becomes : 



Z° N {T,H) = V^Z^^ + Ml-^Z^i^H) 

\ 2(JV-1) _ -v 2(/V-l) 

+4(1 - r)r- 3 cosh 2 (/?//) x + 2 — 2 (15) 

A_|_ — A_ 

Similarly the expression (7) for the open chain partition function with odd number of 
spins takes the form: 



Z° 2N+1 (T, H) = v^Vi + 2(1 - r)cosh (/3H) 

\ 2N > 2JV 
A + - A+ 

Exact expressions for the thermodynamic functions can be calculated by well known 
methods [3]. 

Ising model on Fibonacci chain 

A Fibonacci chain can be inflated by two bonds L(large) and S(small) by the inflation 
ruleL — > LS, S — > L.The chain can be represented by the sequence: 

L — > LS — ► LSL — ► LSLLS — ► LSLLSLSL — ► (17) 

In this case the interaction strengths in the Hamiltonian (1) e i i+1 = e for long bonds 
and e^j+i = e for the short ones where the bonds are arranged according to the Fibonacci 
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sequence (17). The corresponding partition function of the nth generation Fibonacci chain 
having N spins with N — 1 bonds is given by: 

Z° N (F) = TrPPPPP P (18) 

where for long bonds the transfer matrix P is given by eqn.(4) and for short bonds the 
transfer matrix P is given by eqn.(4) with r replaced by f = r\ r Henceforth Z^(F) and 
Zx(I) will represent partition functions for Ising models on an open Fibonacci chain and 
on an open regular lattice respectively. The expressions for the partition functions with 
odd and even number of bonds take the same forms as shown in eqns.(6) and (7) with x^s 
given in eqn.(8) for long bonds whereas for short bonds we replace r by f in eqn.(8).The 
explicit expressions for the partition functions for open and closed chains are: 

Z° N {F) = r—f-Tr{\ + x x ){\ + x 2 ) (1 + x 2N _ 1 )(l + A ) (19) 

Z 2N+1 (F) = r-f*Tr(l + Xl )(l + x 2 ) (1 + x 2N )(l + X)a, (20) 

Similarly for closed chain eqs. (18), (19) and (20) take the following forms: 

Z C N (F) = TrPPPPP P (21) 

Nt JVq 

Z C 2N {F) = r—r*Tr(l + x 1 )(l + x 2 ) (1 + x 2JV _ 1 )(l + x 2N ) (22) 

Nt Nq 

Z C 2 N+1 (F) = r-f*Tr{\ + x x )(l + x 2 ) (1 + x 2N )(l + x^)^ (23) 

where Nl,N$ are number of long and short bonds in a particular sequence. 
Now P is related to P through the following equation: 

P = Vf (1 - ^)<Ti + (24) 

Using eqn.(24) in eqns.(19) and (20) the Fibonacci partition function for any genera- 
tion can be written in terms of open Ising partition functions as follows: 

N 

Z° 2N (F) = h (e, c) + J2h2i(e, e)Z° 2i (I) (25) 
i=i 

N 

Z° 2N -,{F) = k(e, e) + e)Z° 2i _ x {I) (26) 

i=i 

where Z 2i (I) and Z 2i _ x {l) are given by eqns.(15) and (16) respectively. We observe that 
the quasiperiodic nature of the Fibonacci chain is encoded in the functions h(e, e) and 
Z(e, e). Though for small generations these functions can be derived exactly still we could 
not find out their general forms. 

To circumvent this difficulty we study the recurrence relations among the partition 
functions of different Fibonacci generations. A servey of different Fibonacci generations 
depicted by eqn.(17) shows a symmetric pattern in terms of the number of bonds ,viz., 
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Pi : P(odd) 
P 2 : PP{even) 
P 3 : PPP(odd) 
P 4 : PPPPP(odd) 
P 5 : PPPPPPPP(even) 
P 6 : PPPPPPPPPPPPP(odd) (27) 

and so on. 



Recurrence relation among partition functions 

Let P n -2 be the (n — 2)th Fibonacci generation with even number of bonds. This 
automatically ensures that the previous as well as the next two consecutive generations 
will have odd number of bonds . The recurrence relation for Fibonacci generations is 
given by: 

Pn = Pn-lPn-2 (28) 

Now adding a term D n _ 2 P n _3 in the above equation gives 

Pn + P J n-2Pn-3 ~ Pi-lPi-2 + -Dn-2-Pi-3 (29) 

where D n _ 2 = Pet(P„_ 2 ).The following operations are applied sequentially on eqn.(29): 

Substitute P n _ 2 ~ x P n _\ in place of P„_3 on the right hand side and finally use Cayley- 
Hamilton theorem to get the usual trace map relation on the Fibonacci lattice: 

TrP n = TrPn^TrPn-2 - D n . 2 TrP n - 3 (30) 

The above equation will be necessary for calculating recurrence relation among dif- 
ferent Fibonacci generations. For this purpose we must understand symmetry properties 
of Fibonacci chain. Inspecting different generations of the Fibonacci chain it reveals that 
if the total number of bonds N of a particular generation is odd then there is a mirror 

reflection symmetry arround the (^Y^Jth bond ;except the last two bonds. If the special 

bond arround which mirror symmetry occurs is a short (long) one the Fibonacci generation 
will have equal number of odd and even short (long) bonds . However if the total number 
of bonds iV is even,the mirror reflection symmetry is arround a cluster of two successive 

long bonds at the (jf^th and (^j^jth positions of the chain . So "Mirror reflection 

symmetry" is a characteristic property of a Fibonacci chain. 

The nth and (n ± 3)th generations have the mirror reflection symmetry property 
arround the same kind of bond with last two bonds interchanged , while the (n ± 6)th 
generations are topologically same as the nth one. 

Using recurrence relation (28) we can write 

Pn-2pi-3 — Dn-2P n -2 1 P n -1 (31) 

Using Cayley-Hamilton theorem on the right hand side of eq.(31) we get: 

Dn-2Pn-3 — \T^Pn-2)Pn-l — Pi-2pi-l (32) 

Multiplying eq.(32) by P from the right and taking trace we obtain: 
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D n _ 2 Z c n _ 3 = (TrPn^Z^ - Tr{P n . 2 P n ^P) (33) 
In a similar fashion we obtain : 

D n _ 2 Z° n _ 3 = ^ p n-i)K-i ~ rr(P n _ 2 P„_ 1 P ) (34) 

The expression P n _ 2 P n _i in eqns.(33) and (34) is similar to P n = P n -iP n -2 with 
last two bonds interchanged, i.e., both of them have the same mirror symmetric part 
tt n . Therefore eqns.(33) and (34) can be written as: 

D n _ 2 Z c n _ 3 = Z c n _ x (TrP n ^ 2 ) - Tr(Q n PPP) (35) 

D n _ 2 Z°_ 3 = Z° n _ x (TrP n ^ 2 ) - Tr(n n PPP ) (36) 

The transfer matrix has the property that P = P T and P = P T .If such transfer 
matrices are arranged in a mirror symmetric fashion then the resulting matrix (Q n ) will 
have the following properties: 

i) Off diagonal elements are same i.e.,(u; n )i 2 = (uj n ) 2 i 

ii) Diagonal elements are not same but satisfy the condition: 
(u n )u(p,q) = (v n )22(q,p); where p = e pH ,q = e~ m . 

Thus the matrix Q n in eqns. (34) and (35) is of the form: 

q _ ( (^n)ll (^n)l2 \ 

Eqns. (35) and (36) can be written explictly in the following way: 



D n _ 2 Z c n _ 3 = Z^ 1 (TrP n _ 2 )-rv^[(a; n )ii(y + y) + K)i2^ + ^ + ^^) 
+K) 22 (s+y)] (37) 



D n _ 2 Z°_ 3 = Z° n ^ x (JrP n _ 2 ) - Vrr [(co n ) u (y +pu) + (u n ) 12 (u + v + px + qy) 
+(u n ) 22 (x + qv)\ (38) 

where we have used 



PP = yfif 

with 



u y 

x v 



x = 2 + 2 y = p + i u = 1 + V-_ an d v = 1 + a = . 
Eliminating TrP n _ 2 from eqns. (37) and (38) we have: 

yV n -i+puV' n -i xV n _ 1 + qvV'n-! 

-(W„)n + ; r-- — ; -7- (W„) 22 + (U> n ) 12 



(u + v)V n -i + (px + qy)V' n -i (u + v)V n -i + (px + qy)V' n -i 

r-v/F (m + v)K-i + (pa; + ^)V / „-i 
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where 

K-i = ZJ ^r~ K-i and n-i - - — . 

To solve for (co n )ii,(u> n )i 2 and (w n ) 2 2 we need another two equations. These equations 
are obtained from the usual formulae: 

Z<; t = Tr(n n PPP), Z° n = Tr(n n PPP ). 

The explicit forms of these two relations are: 

x + r — y + — 

' (^n)ll + r ax+m {u n )22 + (^n)l2 



r r 



Z n 



r\fru + v + 



qx+py 



r 



x + pu y + gv 

W)n H ; ; ; (^22 + (wJi2 



■u + v + gx + p?/ it + v + gx + 

^° 1 



y/rr u + v + qx + py 

By elementary calculation one obtains: 



r — v n 

I \ / \ x n— 1 x n—l J - y n—2 a 

K)n(z,y) = p 77 77; t x — 7= x A n _! 

1 n _iA' n _i - 1 n-iA n _i rvr 

1 1 7C 



and 



Wu(l '» )= (^ + W)(l-i) 
where 



1 = x V n - (1 - -) ( pu(a;„)ii + qv(u n ) 2 2 



(40) 



(41) 



x {Y n _ x Z° n K' - V'^K) (42) 



v^rViAVi-rViAn-i v n n 



/ \ / \ A ra _i — A' n _i -Dn-2 A 

K) 22 (z,y) = -7 — T7 77; — 7 — X — X A n _! 

l„_iAVi-l ti-iAn-i rvr 

+ 4=77— 1 r/ 7 x iV n ^Z° n K' - Y'^K) (43) 

Vrr r n _iA'„_i - P n _iA n _i 



(44) 



a ( yVn-l+puV'n-! X + f 

An ~ l[X ' y) {u + v)V n ^ + {px + qy)V' n ^ u + v + ^ [ ^ } 
n ~ lMJ (« + v)K_i + (pa; + gy)y' n _i u + v + gx + py 1 j 
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T" n -i(x,y) 



xV n -i + qvV'n-x 



y + qv 



(u + v)V n - X + (px + qy)V' n _i u + v + qx + py 

ryo rye <JC <JO 



A„_i(a;,y) 



K(x,y) 



(u + + (px + gy)y n _i 

1 



u + v + 



qx+py 



K'(x,y) = 

u + v + qx + py 

Eliminating D n _ 2 from eqns.(37) and (38) we obtain: 



(48) 
(49) 
(50) 
(51) 



TrP, 



r^r~ x \/ n -3 



n-2 



70 rye ry c ryo 



O-xy + Pxy -^J- 3 ) (^n)ll(x, y) 



K-3 
V n -3 



( 1 1 * r 

+ xj + 



\^n)22 



x,y) 



rr(l - ^)(<?x + py) 



« + f ) + (px + gj/) ^, n 3 ) 1 



v; 



n— 3 



where 

a. 



a 



y - 

x — 



u+v 
py+qx 

u+v 



X pu, 
x gt>, 



Pxy 



(x-y)(q-p) 
py+qx 

_ (x-y)(q-p) 



x pu 
x gi> 



(52) 



X V py+qx ^ " ' xy py+qx 

In general nth and (n±2)th generations have same arrangement of the last two bonds 
appart from their respective mirror symmetric parts. That is why TrP n and TrP n ±2 will 
have similar expressions. Since we have assumed P n = Q n PP it follows from the expression 
(27) that P n _3 = Q n -sPP. Therefore proceeding in a similar way as above we get: 



TrP 



r\Jr x V„- 



n— 4 



n-3 



70 rye 



J n-2^n-4 



Qyx Pyx T/ - : (^n-l)ll(yj x ) 



V r 



n— 4 



+ (a'yx + P'yjy^ ) (wn-i)n(y, x) 



n-4 



^(l - i)(px + qry) 



w + f ) + (qx + py) 



V 



n—4 



K-4 



)K-1 



(53) 



Using eqns. (52) and (53) and similar expressions for TrP n , TrP n _\ in the trace 
map relation (30) we obtain the following recurrence relation among partition functions 
of different Fibonacci generations as: 



K-i 



ZO rye ry c ryo 

n+l^n-l ~ ^n+l^n-l 



Otxy ~\~ ftxy~ 



V n -^J (uj n+2 )u(x, y) + (oL xy + (3' xy ^-^J {uj n +2)22{x, y) 



^rr(l - l)(qx+py) 



u + v + (px + qy) ■ 



Vn-l 



n+2 



r\fr x V r 



n-2 



^n^n-2 ~ ZnZn-2 
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V, 



x ■ 



n-3 



Zo yc yc yo 

n-l^n-3 ~ ^n-l^n-3 

1 

rf(l - i)(px + qy) 

V'n-3 



( ttj/a; Pyx T , 



D n -2V n -4 



Vn-3 



u + v + (qx + £>2/) 



W 2 2(x,y) 



n-2 



V n -2 



v r 



n+1 



X 



n— 3 



70 yc 
z 'n-2 z 'n-4 



z 'n-2 z 'n-4 



(a yx + Pyx-^-^j (wn-i)n(j/, x) + (a' ya: + 



V'n-3 

m + v + (px + qy) 

V n -4 



{u n ) u (x,y) 

V ' n -3 



K-3 

122(1/, x) 



rr(l - \){px + gy) 



m + v + (qx + py) 



V n - 



(54) 



The above equation reveals the reccurence relation among the partition functions of 
different Fibonacci generations from (n — 4)th to (n + 2)th.The partition functions have 
entered in the above equation through the quantities given by equations from (42) to (51) 
.This is in conformity with the symmetry properties of the Fibonacci chain. 



Conclusion 

We have found an exact expression for the partition function of an Ising model on a 
reguler lattice in presence of magnetic field with open boundary conditions. This always 
includes closed partition functions because of the fact that out of four different spin 
configurations at the end points of the chain two configurations ft an d II satisfy closed 
boundary conditions. We have also shown that "Mirror Symmetry" is a characteristic 
property of all Fibonacci generations . This is the property which leads to a recurrence 
relation among the partition functions of different Fibonacci generations from nth to 
(n ± 6)th . Since nth and (n ± 6)th generations are topologically same , one must go 
through six times decimation renormalization group procedure to find scaling forms of 
thermodynamic functions. 
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